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ABSTRACT 

We constrain the deviation of adiabatic evolution of the Universe using the data on the Cosmic 
Microwave Background (CMB) temperature anisotropies measured by the Planck satellite and a sam¬ 
ple of 481 X-ray selected clusters with spectroscopically measured redshifts. To avoid antenna beam 
effects, we bring all the maps to the same resolution. We use a CMB template to subtract the cosmo¬ 
logical signal while preserving the Thermal Sunyaev-Zeldovich (TSZ) anisotropies; next, we remove 
galactic foreground emissions around each cluster and we mask out all known point sources. If the 
CMB black-body temperature scales with redshift as T{z) = Tq{1 + we constrain deviations 

of adiabatic evolution to be a = —0.007 ± 0.013, consistent with the temperature-redshift relation of 
the standard cosmological model. This result could suffer from a potential bias 6a associated with 
the CMB template, that we quantify it to be |(5a| < 0.02 and with the same sign than the measured 
value of a, but is free from those biases associated with using TSZ selected clusters; it represents the 
best constraint to date of the temperature-redshift relation of the Big-Bang model using only CMB 
data, confirming previous results. 

Subject headings: Cosmic Microwave Background. Cosmology: theory. Cosmology: observations 


1. INTRODUCTION 

Adiabatic expansion and photon number conservati on have produ ced a Cosmic Microwave Background (CMB) with 
a black-body temperature of Tq = 2.725 ± 0.002K (jFixsenI 120091 1 that evolves with r edshift z as T(z) = 1^ (1 -I- 
z). This temperature-r edshift relation is an important test of the Big-Bang paradigm (jAvgoustidis et a1l2012tl and 
of spatial homogeneity (jClarksonll2012tl . Model s like decaying yacuu i n energy density and OTavitationa.1 ‘adiabatic’ 
photo n creation predict a more general scaling (|Overduin et al.l[19981 : iMatviasek et al.l [19951 : |Pu\^ 12004 iJetzer et al.l 
1201111 . An imprint on the T(z) relation can be pro duced if the period of accelerated expansion is driven by a phase 
transition (|Mortonson et al1l2009l : iNunes et al.ll2009ll . The non-conservation of the photon number density changes the 
temperature-redshift and the distance duality relations. Two functional forms ha ve been considered in the literature 
T{z) = To(1-|-2:)^““ (jLima et al.ll2000ll and T{z) = To{l + bz) (iLoSecco et al.l[2001^ . with a and b constant parameters; 
at low redshifts, the differences between both functional forms are small and only the first scaling is usually tested. 
In most scena rios, deviations of adiabatic evolution are associated with distortion s of the CMB black-body spectrum 
l|Chlubal[201411 and are strongly constrained by the current FIRAS upper limit of iFixseiJ (|2009il . A departure of the 
standard temperature-redshift relation would represent an important challenge to the current cosmological model and 
it represents a test of these alternative scenarios. 

The earliest measurements on CMB black-body temperature evolution were o btained using the relative population s 
of atomic fine-structure levels which are excited by the background radiation (|Songailalll994 l^rianand et al.ll200nll . 
iNoterdaeme et all (|2010l l measured T(z = 2.69) = 10.5^® gK using the rotational excitation of CO molecules in quasar 
spectral lines, value compatible with 10.06K, the black-body temperature expected at that redshift with the standard 
temperature-redshift relation. Using quasar spect ral lines, the best constraint on deviations from adiabatic evolution 
at present is a = 0.009 ± 0.019 (jMuller et al.ll2013h : they also obtained a very stringent individual measurement of the 
CMB temperature: T(z = 0 .89) = 5.08±0.10K. Recently the Planck Collaboration obtained a much tighter constraint 
([Planck Collaboration! |2M^ by including data on large scale structure to the CMB data but their measurement 
a = (0.2 ± 1.1) X 10“^ applies to models where the deviation from adiabatic evolution starts at the last scattering 
surface. The temperature-redshift relation has also been probed using the Thermal Sunyaev-Zeldovich effect (Sunyaev 
& Zeldovich, 1970; hereafter TSZ). The TSZ anisotropy induced by clusters of galaxies along the line of sight h is 
AT(h) = TqG{x)Yc. The Comptonization parameter is defined as Yc = {kB/iTiec'^) J T^dr, where dr = arUedl is 
the cluster optical depth, ne,Tf, are the electron density and temperature evaluated along the line of sight Z, ax the 
Thomson cross section, ks the Boltzmann constant, me the electron mass and c the speed of light. This effect depends 
on the frequency of observation i/q as x = hi'Q{l + z)/kBT, where h is the Planck constant and T = T{z) the CMB 
black-body temperature at the cluster location. In the non-relativistic limit, G{x) = a;coth(a;/2) — 4. This spectral 
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dependence is different from that of all known foregrounds, making the TSZ effect an effective tool for detecting clusters 
as well as a potential probe of the redshift evolution of the background temperature. 

If the Univ erse evolves adiaba tically then x is independent of redshift. We will test deviations fro m adiabatic ev o lution 
assuming the lLima et al.l (1200(11 1 parametrization, that is, x = hvQ{l + z)°'/U bTq and G = G{v,a). iFabbri et al.l (jl978[ l 
proposed to constrain a by measuring the zero cross frequency of clusters at different redshifts that, for adiabatic 
evolution, occurs at v ^ 217GHz. The measurement of the cross-over frequency is problematic since the TSZ is 
inherently weak and could be dominated by uncertain systematics. As an alternative, early studies fit the T SZ signal 
at diff erent frequencies to measure the function G'(:/, a). We shall denote this procedure the fit method. iReohaelil 
(jl980ll suggested the ratio method, constraining a by using the ratio of the TSZ anisotropy at different frequencies, 
R{vi,V 2 ,oi) = G(vi, a)/G{ v 2 , a). By taking ratios, the dependence on the Comptonization parameter is removed and 
there is no need to account for model uncertainties on the gas density and temperature profiles. On the negative side, 
the analysis is mo re complicated since the distribution of the temperature anisotropy ratios is highly non-gaussian 
(jLuzzi et al.ll2009ll and the Kinematic Sunyaev-Zeldovich (Sunyaev & Zeldovich, 1972; hereafter KSZ) generated by 
cluster peculiar velocities is a component always present at the cluster location that can not be separated from the 
cosmological signal. For small number of clusters the effect needs to be taken into account while it is expected to be 
less important for large cluster samples. iBattistelli et al.l (l2002h found a = — from the TSZ measurement of 
the Coma and A2163 clusters. Later, from the data on just 13 clusters, iLuzzi et ahl (12009( 1 set up an upper limit of 
a < 0.079 in the range z = 0.023 — 0.546 at the 68% confidence level. 

In 2013 the Planck Collaboration released their first maps of the temperature fluctuations of the CMB sky using 
15.5 months of data (iPlanck Collaborationll2014ajl . The Low Frequency Instrument (LFI) produced three maps with 
frequencies 30, 44 and 70 GHz and angular resolution from 33 to 13'. The High Frequency Instrument (HFI) was 
sensitive to a wider range of frequencies, from 100 to 857 GHz and scanned the CMB sky at much higher angular 
resolution, from 10 to 5'. Due to its large frequency coverage, high resolution and low noise, the Planck satellite is an 
optimal instrument for detecting the TSZ distortion of the CMB spectrum induced by clusters. To clean the CMB data 
from foreground contributions, the Planck Collaboration also released maps of (a) thermal dust and residual Cosmic 
Infrared Background emission, (b) synchrotron, free-free and spinning dust emission, ( c) CO contributions, that ar e 
most important at 100, 217 and 353 GHz and (d) maps of dust temperature and opacity (|Planck Collaborationll20143l . 
Additional products required to analyze the data suc h as masks and nois e maps were also released. Using Planck data 
and a sample of 813 clusters up to redshift z ~ 1, iHurier et al.l (|2014H obtained a = 0.009 ± 0.0 17. Their s ample 
included X-ray and TSZ selected clusters. Using clusters detected with the South Pole Telescope, iSaro et abl (|2014(1 
measured a = 0.017to'o28j also compatible with adiabatic evolution. 

In this article we will apply the techniques developed in ide Martino et ^ (1201211 to test the standard temperature- 
redshift relation. Our cluster catalog comprises 782 clusters selected from the X-ray ROSAT data and with well 
measured properties; of those, 481 will be used in this study. The outline of this paper is as follows: In Sec. 2 we 
describe our data and pipeline; in Sec. 3 we present our results and in Sec. 4 we summarize our main conclusions. 

2. DATA AND DATA PROCESSING 

We constrain the temperature-redshift evolution of the background temperature using a cluster catalog selected 
from ROSAT data and the Planck Nominal maps released in 201 30. We will only use HFI data, originally released in 
Healpix format with resolution Ngide = 2048 (iGorski et al.ll20(i5ll . 


2.1. Cluster Sample. 

Our catalog contains 782 clusters, of which 623 are outside the Planck PCCS-SZ-Union mask (jPlanck CollaborationI 
(see below). They have been selected from three X-ray flux limited surveys: the ROSAT-ESO Flux Limited 
X-ray catalog (REFLEX, Bohringer et al. 2004), the extended Brightest Gluster Sample (eBGS, Ebeling et al. 1998 
and 200 2) and the Glusters in the Z one of Avoidance (GIZA, Ebeling et al. 2002). Details of the combined catalog are 
given in iKocevski fc Ebelind (I2006D . All clusters have well measured positions, X-ray fluxes and luminosities in the 
[0.1,2.4]KeV ROSAT band, spectroscopic redshifts and angular e xtents of the X-ray emitting region. The catalog also 
lists the X-ray temperature derived from the Lx — Tx relation of lWhite et al.l (|1997r i and the core radii Tc and central 
electron density Uep obtained by fitting a ,8 = 2/3 model to the ROSAT data. We compute the rad ius at which the 
mean overdensity of the cluster is 500 times the critical density, rsoo, from the r^oo ~ Lx relation of iBdhringer et ^ 
(|2007ll . We define M^oo as the mass enclosed within a sphere of radius 7-500 and the angular size ^500 = 7 - 500 /d^, where 
dA{z) is the diameter angular distance of each cluster in the ACDM model. In Fig. [T]we represent the redshift and 
mass distribution of the clusters in our catalog. 

The clusters in our sample have masses in the interval M 500 = [0.2 — 10] x l O^^Mfn and are loc ated at redshifts 
z < 0.3, relatively lower compared with other catalogs. The catalog used byjHurieret^ (120141 1 contains X-ray 
and Planck selected clusters. Compared with our cata log, the MCXC sampl e (iPiffarett i 2011 ) includes clu sters from 
sources that we ha ve not considered such as NORAS (|B5hringer et al.l[200011 . SGP ([Cruddace et al.ll200§l and NEP 
([Henry et al.ll200(ill . Also, our catalog does not include clusters from X-ray pointed surveys which tend to find mainly 
low-mass systems but out to high redshift. For redshifts z < 0.3 and masses M 500 > all our clusters except 

ten are listed in the MCXC sample. In the same mass and redshift range the Planck SZ catalog contains 555 clusters. 


® Planck data can be downloaded from http://www.cosmos.6sa.int/web/planck 
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Of those, only 239 are at less than 10' away from members of our sample. An important diffe rence with the latter 
catalog is that the Planck Collaboration e stimates M^nn us i ng the scaling relation M 500 —PA.soo of lArnaud et al.l (|2010(1 
while we u s e the r^nn — Lx relation from iBohrineer et al.l (|2007f) and the resulting masses differ, on average, by 30%. 
iSaro et id] ()2014l l used a sample of 158 clusters with redshifts z < 1.35, selected by their TSZ signal. When cluster 
candidates are identified using CMB data the adiabatic evolution is assumed, biasing cluster selection towards those 
candidates that mimic this behavior. The South Pole group verified that this effect was not significant. Nevertheless, 
due to their selection criteria, negative temperature fluctuation s were more like ly negative than positive and taking 
into account this second effect increased their error bars a 30% (ISaro et al.ll2014ll . Therefore, it is important to verify 
the results of other groups using only X-ray selected clusters since they will not be affected by these biases. In Table |T] 
we present the mean redshift, angular scale. X-ray luminosity and mass averaged over the full cluster sample and in 
several subsamples selected acc o rding to X-ray luminosity mass and redshift. For better comparison we chose the 
redshift bins as in iHurier et al.l (j2014l l: within each redshift bin, we select all clusters with M 500 > 2 x IO^^Mq or 
Lx > 2.5 X lO^'^ erg/s. In the table we do not quote a redshift bin in the mass and X-ray luminosity subsamples when 
it coincides with the bin of the full sample. 


Subset 

Ncl 

z 

^500 

(arcmin) 

Lx 

(10‘^^erg/s) 

M500 
(1014 Mq) 

All Clusters 

481 

0.106 

12.3 

2.35 

3.1 

0.0 < z < 0.05 

32 

0.035 

21.2 

0.51 

1.2 

0.05 < z < 0.10 

186 

0.074 

12.5 

1.34 

2.4 

0.10 < z < 0.15 

114 

0.123 

8.82 

2.20 

3.3 

0.15 < z < 0.20 

83 

0.169 

7.64 

4.17 

4.8 

0.20 < z < 0.25 

46 

0.222 

6.58 

6.13 

5.8 

0.25 < z < 0.30 

20 

0.274 

6.15 

10.1 

7.6 

M500 > 2 X IO^^Mq 

397 

0.134 

10.49 

3.40 

4.1 

0.0 < z < 0.05 

20 

0.039 

26.5 

1.67 

3.1 

0.05 < z < 0.10 

121 

0.078 

13.4 

2.01 

3.2 

0.10 < z < 0.15 

107 

0.124 

8.85 

2.28 

3.4 

Lx > 2.5 X lO^^erg/s 

201 

0.152 

13.51 

5.32 

5.7 

0.0 < z < 0.05 

3 

0.037 

38.4 

3.47 

5.2 

0.05 < z < 0.10 

25 

0.077 

16.6 

4.28 

5.6 

0.10 < z < 0.15 

36 

0.129 

9.59 

3.51 

4.6 

0.15 < z < 0.20 

71 

0.171 

7.72 

4.51 

5.1 


TABLE 1 

Average properties of different cluster subsamples from the general catalog. Clusters have been selected by luminosity 

AND MASS WITHIN EACH REDSHIFT BIN. 


2.2. Foreground Cleaned Planek Nominal Maps. 

The nine frequency maps released by the Planck Collaboration in 2013 contained foreground emissions in addition to 
the intrinsic CMB temperature anisotropies and the instrumental noise. Since the different frequencies have different 
angular resolutions, we bring all maps to the common resolution of 10', the lowest of the HFI channels. Due to 
their lower resolution and higher instrumental noise we will not use the LFI data. We r emove the cosmological and 
KSZ anisotropies by subtracting the LGMCA CMB map from iBobin et al.l (|2013L 1201411 smoothed to the same 10' 
resolution. This map was constructed from the latest Wilkinson Mierowave Anisotropy Probe 9yr and first Planck 2013 
data releases. The data were combined using a component separation method based on the sparsity of the foregrounds 
in the wavelet domain that requires frequency information to remove the TSZ effect. The resulting map has low dust 
contamination at £ < 1000 and the KSZ anisotropy is preserved while the TSZ signal is removed. In this respect, 
it is more useful than the Planck reconstructions of the CMB anisotropies, like SMICA, that contain substantial 
TSZ residuals. Next, we clean the foreground emissio n around each cluster. At fr equencies higher than 100 GHz, 
thermal dust emission dominates over most of the sky (|Planck Collaborationll^lddl . This contribution is commonly 
described as a modified blac k-body spectrum with a power-law emissivity e,^ oc with a slope « 1.5 — 1.8 
(jPlanck Collaborationll2014bll . independent of frequency but varying across the sky even on scales as small as 5'. The 
dust contribution is largest at 857 GHz and dominates over the cosmological signal. Then, we can use this channel as a 
template for thermal dust to clean patches P(;/, x) centered on each cluster position x and on each maps of frequency 
v. When the dust emission parameters and frequency dependence are known, one can generate dust templates at 
different frequenci es to be subtracted from the data. In this paper we will use a different approach, proposed by 
iDiego et al.l (l2002f) . that minimizes the contribution to the map at frequency v of all the foregrounds that correlate 
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with the 857 GHz channel. Specifically, the weight w{v) minimizes the difference [P(;/,x) — w{i')P{857 GHz,x)]^. 
Weights are computed on a ring around the cluster region, defined as 77. = [9ci , Opatch], and are given by 

., Sxe7^P(Gx)P(857 GHz,x) 

= S,,„|P(857 GH.,x)P ^ <'> 

To avoid overlapping with the TSZ emission from the clusters in our catalog, we masked a region of S^soo around 
them. Excising these pixels changes the weight by less than 2%. We assume that the dust temperature and spectral 
index are constant within the ring. This approximation is more accurate for the less extended clusters. Ring sizes vary 
depending on the cluster extent: 9ci = 26 * 500 , 9patch = 305 oo for all clusters with 6*500 > 20 ', otherwise 9ci = 39^00 and 
9patch = [5 — 7] 6 * 500 . If we would use a ring of the same size for all clusters, the patch would be too large for the most 
extended clusters and the approximation of constant temperature and constant spectral index will not held. On the 
contrary, if 9500 is too small, then the patch would contain few pixels and the statistical weight would be unreliable. 
We checked that the above angular extents yield the smallest foreground residuals around clusters. The procedure is 
robust since weights did not change if 9patch and 9^1 varied by a factor of two. Then, even if the smallest/largest annuli 
could contain a small TSZ contribution from the cluster or nearby clusters, the weights were not affected. 

The first three CO rotational transition lines J=f^-0, J = 2—>-1, and J = 3 —>■ 2 at 115, 231 and 346 GHz, 
respectively, present the largest transmission coefficients making them a significant foreground component in the 
Planck intensity maps. We corrected this contribution using the esti mated CO emission maps pro vided by the Planck 
Collaboration. Three type of CO maps has been made available (iPlanck Collaborationll^lddf) : Type 1 maps are 
available for the three frequencies but are too noisy and are sensitive only to the brightest regions on the galactic plane 
to be useful. Type 2 maps are less noisy but they are only available at 100 and 217 GHz. We corrected these two 
frequencies and found the correction to affect very little the final cleaned patches. Since the transition at 353 GHz is 
the weakest of the three, we expect that the CO contamination at 353 GHz will neither produce a significant effect on 
the final results. 

To reduce the contamination from point sou rces and foreground residuals near the galactic plane, we used the PCCS- 
SZ-Union mask (|Planck Collaborationll^ldd ill. This mask was constructed using the Planck Catalogue of Compact 
Sources (PCCS). It is the union of six masks, one for each HFI channel, and the mask of the Galactic Plane and the 
Magellanic Clouds. Our cleaning method does not give satisfactory results for the faintest clusters, so we restricted our 
analysis to those clusters with a X-ray luminosity in the ROSAT [0.1-2.4]KeV band of Lx > 0.5 x 10"^'*erg/s, reducing 
the total number of clusters to 481. 

To illustrate how effectively our pipeline removes foregrounds and the differences in the final results between low 
and extended and high and compact redshift clusters, in Fig. [5] we show the temperature fluctuation on patches of 
angular size 75' x 75' for two clusters: Coma, with an angular extent of ^500 = 48.1' and redshift z = 0.023, and PSZl 
G355.07-1-46.20 with 6*500 = 9.2' and z = 0.21. The first and third rows correspond to the temperature anisotropies on 
the original Planck Nominal maps at frequencies 100 — 545 GHz. The second and fourth rows show the same regions on 
the foreground cleaned maps. In Coma and PSZl G355.07-1-46.20 we fixed the temperature range to be [—300,300]/rK 
and [—200,200]/j,K, respectively. In the Planck nominal maps, the two circles have a angular radius of 9ci (inner 
circle) and 9patch (outer circle). In the foreground clean maps the radius corresponds to 6500 - The TSZ temperature 
anisotropy is negative at 100 and 143GHz; at 217 GHz, it is greatly reduced and it changes sign at higher frequencies. 
The data at 545 GHz is dominated by foreground residuals and our pipeline produces reliable estimates of the TSZ 
signal only for the brightest clusters. At this frequency, foreground residuals dominate over cluster anisotropies and 
we could not use it to determine/constrain deviations from adiabatic evolution. 

2.3. Testing our Cleaning Procedure 

We checked that the LGMCA map does not contain significant TSZ residuals by computing the mean temperature 
anisotropy at the position of our clusters on discs of different sizes, in units of 05 oo- In Fig. [3^, blue squares and red 
solid circles represent the mean at the actual position of our 481 clusters and the mean at 481 random position in 
the sky, respectively. The error bars are the error on the mean estimated from 100 simulations. For a better view, 
the results are displayed slightly shifted. The averages at the cluster positions are marginally biased towards negative 
values compared with a random distribution. For apertures > ^ 500 , this bias is ~ — l^K or smaller and well within 
the error bar. Since the TSZ is removed from the LGMCA map using the standard (a = 0) frequency dependence of 
the effect, this offset could be systematic and not the mean of some random residuals left at the cluster locations by 
the component separation process. If a 76 0, the TSZ residuals that would remain would change with frequency. If the 
offset were systematic the TSZ effect would be shifted upwards by lp,K after subtracting the LGMCA map, modifying 
the dependence with a but not removing it. We will discuss this point further in Sec. 3.3. 

To illustrate that the TSZ anisotropy has been preserved in our foreground cleaned and CMB subtracted maps, in 
Fig. [ 3)3 we plot the mean temperature anisotropy averaged on a disc of radius 6 * 500 . Open triangles correspond to the 
profile of the Coma cluster, while solid circles to PSZl G355.07-1-46.20. The error bars were computed by placing discs 
with the same angular extent at one thousand random positions outside the clusters in our sample; they are much 
smaller than those of Fig [3^ since only the instrumental noise and foreground residuals contribute. The difference in 
the errors between both clusters are due to their different angular extent. Third, we verified that there are no TSZ 
anisotropy on the rings selected to compute the weights of eq. (P) that could potentially bias the result. Fig. |3t shows 
that the mean temperature anisotropy on the rings around of the two clusters is ~ 10 “^/rK, negligible compared with 
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the errors on the cluster profiles; around other clusters the mean was always < also negligible. To conclude, 

subtracting the LGMCA map from Planck Nominal maps effectively removes the cosmological CMB and KSZ signals 
while preserving the TSZ anisotropy and its dependence with a. 


2.4. Error Bar Estimation 


For each cluster configuration, error bars were computed by evaluating the mean temperature fluctuation on 1,000 
random positions in foreground cleaned maps on a disc with the same angular extent than the clusters in any given 
subsample. The random positions were chosen to be at least 2 degrees away from the location of our clusters to 
guarantee that the random discs will never overlap with them. We removed the foreground contamination using the 
same procedure than at cluster locations. We verified that the mean of the simulations was compatible with zero. The 
correlation matrix between different frequencies was computed by averaging over all simulations 


C{ui,uj) 


{[STivj) - ^l{vi)\[5T{vJ) - 
a{v,)a{vj) 


( 2 ) 


where /.t(r'i) = {5T{vi)), and a{vi) = {\5T{vi) — To compare with previous work, we repeated the process 

and computed the correlation matrix between foreground clea ned maps before su btracting the LGMCA map using 
our full catalog and compared it with the correlation matrix of iHurier et al.l (I2014li . obtained with a different cluster 
catalog. While the latter authors used a fix aperture of 20' for all clusters we chose apertures of size 05 oo for each cluster. 
We found that our off-diagonal terms were slightly smaller; the difference between the two correlation matrices differed 
between 2% for the element C(100GHz,143GHz) to less than 10% for C(100GHz,353GHz), allowing us to conclude 
that our cleaning procedure was comparatively as effective as theirs. 


3. RESULTS 

To estimate/constrain the temperature-redshift relation we used the ratio and fit methods described in the intro¬ 
duction. We derived an X-ray temperature f or all clusters in our catalog from the measured X-ray luminosity using 
the Lx — Tx relation of IWhite et al.l (|1997f) . The cluster temperature allow us to include relativistic corrections 
(|Itoh et al.lll998t iNozawa et al.l 1200611 . If we denote 9e = (ksTe/meC^), the frequency dependence of the TSZ effect 
including relativistic corrections up to fourth order in the cluster temperature, 04 ( 1 '), can be written as 


G 4 {iy) = G{v) 




(3) 


where Yq, ...,V 4 are defined in eqs 2.24-2.29 of INozawa et al.l (|1998f) . The correction is relatively small; for the most 
massive clusters it amounts to a few percent. Gorrections depending on the unknown cluster peculiar velocities were 
not included but for the values expected in the concordance model, they are also negligible. Hereafter, in order to 
simplify the notation we will drop the subindex; it should be understood that relativistic corrections were included 
except when indicated otherwise. Finally, let us mention that as a prior we considered a = [—1,1], subdivided in 2001 
equally spaced steps. 


3.1. Ratio method 

As indicated in Sec. o all maps were reduced to a common resolution of 10' to eliminate the differences in the 
angular resolution. In this statistic, the ratio of the mean temperature anisotropy on a disc of a fixed angular size 
9cl at two frequencies, ST{lyi)/ST{ 1 / 2 ), does not depend on the cluster pressure profile. In the absence of noise, this 
ratio would be equal to the theoretical value at each redshift R{i/i,V 2 ,a) = G{vi,a)/G{h' 2 ,a), providing a direct 
measurement of a . However, instrumental noise and foreground residuals complicate the analysis. As discussed in 
iLuzzi et all (|2009ll . for each cluster j the probability distribution of the ratios, Vj{R), is 


Vj{R{yi,v2,a)) = 


1 


27r(T,y, ( 7 1/ 


X X exp 


[x-ST{vi)Y , [xR{vi,V 2 ,oi) - ST{v 2 )f 


2c^^i 




dx, 


(4) 


where ai, is the error on mean temperature at the cluster location j at frequency v. The likelihood function is £ oc LijPj ; 
we did not include the correlation between the different frequencies so our analysis will underestimate the errors. We 
took five ratios between frequencies from 100 GHz to 353 GHz: 100/353, 143/353, 217/353, 100/143, 217/143. We 
chose in the denominator the channel with the smallest instrumental noise and excluded 217 GHz to avoid dividing b y 
zero when a ~ 0. Temperature anisotropies were eyaluated on discs of radii 9ci = ^soo- Following iHurier et al.l (j2014[) . 
we diyided our sample in six redshift bins of equal width Az = 0.05; then, in eq. ([4]) ST{vi) and R{vi,V 2 ,a) are now 
the ayerages oyer all the clusters in the bin. 

In Fig. |4] we plot the likelihood function for different cluster subsets: In (a), we considered the 201 most X-ray 
luminous clusters, with X-ray luminosity Lx > 2.5 x lO'^^erg/s, in (b) the 397 clusters with mass M 500 > 2 x IO^^Mq 
and, finally, in (c) we plot the likelihood of our full cluster sample. For illustration, only the likelihood of three redshift 
bins are plotted. In each panel, the dotted, dash-dotted and dashed lines correspond to the bins [0.05,0.1], [0.15,0.2], 
[0.25,0.3], respectively, and the solid lines correspond to the full sample. The measured values of a are given in 
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Table [5] within each redshift bin we computed a in different subsamples, with clusters selected in X-ray luminosity 
and mass as indicated in Table [T] in order to test the relative contribution of the different cluster subsamples to the 
final error budget. In all redshift bins, the results are always compatible with zero. Notice that our error bars do not 
necessarily scale with the square root of number of clusters since, as indicated in Sec. 2.2, our cleaning procedure, and 
consequently the amplitude of the foreground residuals, depend on the cluster extent. We checked that the results were 
very similar if the temperature averages were taken on discs of radii 9ci = 205 oo. Using the full sample, we constrain 
the deviation from adiabatic invariance to be a = —0.03 ± 0.06. This estimate is a factor of 2 worse than the result 
obtained by the SPT group (|Saro et al.ll2014ll . even with our underestimated error bars. This is to be expected since 
their maps have a resolution of 1.5' and even if their sample is smaller in size their clusters are, on average, at higher 
redshifts where the effect of the non-adiabatic evolution is most noticeable. 


Subset 

Nal 

O-Lx 


Ncl 

^-^500 

°'“M'500 

N,i 

“®500 

‘^“<’500 

All 

201 

0.02 

0.06 

397 

0.02 

0.06 

481 

0.02 

0.06 

0.0 < z < 0.05 

3 

0.13 

1.0 

20 

-0.05 

0.74 

32 

0.03 

0.67 

0.05 < z < 0.10 

25 

0.44 

0.77 

121 

0.15 

0.41 

186 

0.06 

0.43 

0.10 < z < 0.15 

36 

0.95 

1.27 

107 

0.32 

0.34 

114 

0.26 

0.29 

0.15 < z < 0.20 

71 

0.56 

0.51 

83 

0.07 

0.17 

83 

0.07 

0.17 

0.20 < z < 0.25 

46 

0.01 

0.01 

46 

0.01 

0.01 

46 

0.01 

0.01 

0.25 < z < 0.30 

20 

-0.01 

0.08 

20 

-0.01 

0.08 

20 

-0.01 

0.08 


TABLE 2 

Values of a estimated using the ratio method. Temperatures are averaged on discs of radius 6500 over all the clusters in 

THE redshift BIN. 


3.2. Frequency fit method 

As an alternative method, we constrain the adiabatic evolution of the Universe_by fitting the frequency dependence 
of the TSZ anis otropy, AT{h ) = TqYcG{i't a). The Comptonization parameter Yq is the average on a disc of radius 
0500 - We follow iHiirier et aU (j2014l l and take it as a free parameter. In this case, both Yc, and a are fit to the data. 
We took the flat prior Yc = [Oj 300]/rK divided in intervals of ISYc = 0.15p,K. Since all maps have the same resolution, 
Yc is independent of frequency. We compute the likelihood function as 


N,i 


k—1 


where C(vi,Vj), given in eq. ([2]), is 


C{ui,uj) = 


Uc,feG(l/i,a)] G 


) [^Tk{vj 

/l.OOOO 

0.9258 

0.4603 

0.1435\ 

0.9258 

1.0000 

0.4995 

0.3168 

0.4603 

0.4995 

1.0000 

0.3605 

\0.1435 

0.3168 

0.3605 

1 .0000/ 


(5) 


( 6 ) 


where Vi = (100,143, 217, 353)GHz. Notice that by using the LGMCA map to subtract the intrinsic CMB signal we 
ha ve significant l y redu ced the correlation between frequencies and the variance of the foreground clean maps compared 
to iHurier et al.l (j2014l l. 

In Fig. [S]we present the likelihood function, marginalized over the Gomptonization parameter, for the same cluster 
subsamples and redshift bins than in Fig. |4l with lines following the same convention. The numerical results are given 
in Table 131 in this case we find that a = —0.007 ± 0.013, compatible with zero. We verified that changing the prior to 
Yc = [—300, 300]/lK, that would allow for foreground residuals to change the sign of the TSZ effect, did not modify 
the error bars, indicating that the data prefers negative values in the Rayleigh-Jeans and positive values in the Wien 
part of the spectrum. To illustrate the accuracy of our results, in Fig. |6]we plot the average temperature anisotropies 
of all the clusters the redshift bins of Table [Hand their errors, for the different frequencies. The solid line corresponds 
to the best fit YcG{v^ a) without relativistic corrections and the dashed line to the best fit including these corrections. 
Notice that only at z > 0.2 this correction can be distinguished from the non-relativistic effect. We also indicate the 
per degree of freedom of the best fit model. 

Our result represents a 30% improvement over iHurier et al.l (l2014li and it was obtained with a smaller sample of 
clusters located at a lower mean redshift. In Fig. [7] we present a more thorough comparison. The shaded regions 
correspond to the 1 and 3a errors given in Table El Notice that in the same redshift range, the measured values are 
in an excellent agreement and are fully compatible with an adiabatic evolution of the Universe, the main difference 
being that our error bars are significantly smaller. The reduction on the uncertainty is due to using an accurate GMB 
template to subtract the cosmological signal. 

We can use the upper limit on a given above to cons t rain t he phenomenological model of iLima et al.l (|2000D using the 
temperature-redshift relation derived bv iJetzer et al.l (|2011l) . The constrain is given in terms of the effective equation 
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Subset 


CtLx 

(J Ot T 

Nd 

^■^500 


Nd 

“^500 


All 

201 

-0.013 

0.014 

397 

-0.006 

0.013 

481 

-0.007 

0.013 

0.0 < 2 < 0.05 

3 

-0.27 

0.23 

20 

-0.21 

0.17 

32 

-0.16 

0.16 

0.05 < z < 0.10 

25 

0.00 

0.18 

121 

0.03 

0.09 

186 

-0.08 

0.09 

0.10 < z < 0.15 

36 

0.08 

0.20 

107 

-0.01 

0.09 

114 

-0.03 

0.07 

0.15 < z < 0.20 

71 

0.03 

0.05 

83 

0.05 

0.04 

83 

0.05 

0.04 

0.20 < z < 0.25 

46 

0.01 

0.02 

46 

0.01 

0.02 

46 

0.01 

0.02 

0.25 < z < 0.30 

20 

-0.03 

0.02 

20 

-0.03 

0.02 

20 

-0.03 

0.02 


TABLE 3 

Estimated values of a and their uncertainties using the fit method. 


of state of a hypothetical decaying dark energy model. From the bounds on the deviation from adiabatic invariance 
we find Weff = —1.005 ± 0.008, fully compati ble with the conco rdance ACDM model. This constrain represents a 20% 
improvement over the latest measurement of ISaro et al.l (|2014l i. 

3.3. Possible systematic effects 

Our estimates do not take into account the effect of subtracting the TSZ effect from the LGMCA map using the 
G{v, a = 0) frequency dependence. Fig. [3^ shows that the average temperature anisotropy at the cluster locations has 
^ — 1/xK residual compared with the same measurement at random positions in the sky averaged over 100 realizations. 
This difference is well within the errors and compatible with sample variance. We checked that the average on the 
six redshift bins was random, oscillating around the mean, indicating that there was not such systematic effect. 
Nevertheless, if the redshift evolution has a ^ 0 then there would exist a residual TSZ at the cluster location and 
we need to consider the biases introduced if the — 1/iK difference was in fact systematic. Since the LGMCA map is a 
combination of WMAP and Planck data at different frequencies, when subtracting the same LGMCA map to all Planck 
Nominal maps, the residual TSZ anisotropy would shift the zero cross frequency but will not change the dependence 
on a. The magnitude and sign of the temperature shift is difficult to predict since it would depend on the cluster 
redshift, mass and extent and on the LGMCA mixing matrices estimated for a set of input channels on a patch of data 
at a given wavelet scale. We reanalyzed the data subtracting a fix 1/iK from the measured temperature anisotropies of 
all the clusters and at all frequencies and obtained that the value of a decreased by 8a = —0.02 in the fit method and 
remained within the errors in the ratio method; i.e., we would be obtaining a value closer to a = 0 by this amount. If 
Of > 0 then the residual TSZ at the cluster location would be positive and so would be the bias in a. Therefore, by 
using the TSZ frequency dependence of an adiabatically evolving Universe, our method could be masking the effect of 
a non-adiabatic evolution. Taking into account this effect, our final constraint would be a = —0.007 ± 0.013 (—0.02) 
where the parenthesis indicates the systematic contribution. At the 2(T level, this result is both compatible with all 
the upper limits from TSZ and spectral lines. 

4. CONCLUSIONS. 

We have constrained the deviation of the adiabatic evolution of the GMB black-body temperature applying two 
different estimators to a sample of X-ray selected clusters and using foreground cleaned Planck Nominal maps. By 
not including clusters selected b y their TSZ signatu re, we avoid biasing our sample to those clusters that are closest 
to adiabatic evolution. Following iHurier et al.l (1201^ we distributed our cluster in six bins of redshift. The constrains 
using the ratio and fit methods were a = —0.03 ±0.06 and a = —0.007 ±0.013, respectively; the fit method producing 
statistically more significant results than the ratio method since the latter is very sensitive to small denominators. The 
constrains are weakened if we add a hypothetical systematic effect due to the component separation method used to 
construct the LGMGA map and become a = —0.007 ± 0.013 (—0.02) and a = —0.03 ± 0.06 (—0.02), compatible with 
the results given i n the literature a t the 2 ct confidence level. 

Gompared with lSaro et al.l (I2014D the errors in the ratio method are at least a factor of two worse. While our sample 
contains four times as many clusters, they are located at much lower redshift, where the statistic is less sensitive. 
However, our implementation of the fit method provides the best constraints on a to date using only GMB data. For 
instance, as illustrated in Fig[3 by using the LGMCA reconstruction of the intrinsic GMB temper ature anisotropies to 
remove the cosmological contribution from our maps, we have reduced the errors compared with IHurier et al.l (j2ni4[ i 
and obtained the best bound to date. The constrain is also interesting since our cluster sample has the lowest mean 
redshift. We constrain deviations of adiabatic evolution within z < 0.3, well in the p eriod of accelerated expansion. 
Our results are comparable with a similar analysis carried out bv iLuzzi et al.l (1201511 using a catalog of clusters at 
higher redshifts and fully complements their findings. 

Since our clusters have been selected using X-ray data, they are not affected by possible biases associated with 
selecting clusters using the TSZ effect, providing a consistency check to previous results. We continue to expand 
our X- ray cluster sarn ple to probe higher redshifts, where systematic biases would be larger and could be estimated 
better. iGhlubal (120141) has argued that low redshift energy injections that could produce a non-adiabatic evolution of 
the Universe are strongly constrained by the upper limit set by FIRAS. Nevertheless, and independently of theoretical 
expectations, measurements of the temperature-redshift relation of the GMB black-body temperature provide a strong 
consistency check of the current Big-Bang paradigm. 
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Fig. 1.— 


Histograms with the mass (a) and redshift (b) distributions of the 481 


clusters used in this study. 
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Fig. 2.— Original Planck Nominal and foreground cleaned patches centered on the position of the Coma [A1656] and the PSZl 
G355.07+46.20 clusters at the HFI frequencies 100-545 GHz. Blue/red colors correspond to negative/positive temperature fluctuations. For 
Coma, the patch is 7^ x 7*^ in size and the temperature range is [—300, 300]/.iK; for PSZl G355.07+46.20 the size of the patch is 2.5*^ x 2.5*^ 
in the Nominal and 1^ x 1^ in the Cleaned maps and the temperature is within [—200, 200]pK. In the Nominal maps, the inner and outer 
circles have radius Od and Opatch^ respectively; in the Cleaned maps, the angular radius is ^500 for tho given cluster. 
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Fig. 3.— (a) Mean temperature anisotropy and its error evaluated on the LGMCA map on discs of different radius, in units of the ^soo- 
Blue open squares, red solid circles correspond to the average at the real position of the 481 real clusters and at the same number of random 
locations on the sky, respectively. The error bars represent the error on the mean, (b) TSZ anisotropy, measured on a disc of size ^500 on 
our CMB removed, foreground cleaned data as a function of frequency for the two clusters represented in Fig. [2] Red circles have been 
shifted for to facilitate their view, (c) Average temperature anisotropy evaluated in the ring [0ci,0patch] for the same two clusters. The 
error for PSZl G355.07+46.20 is too small to be seen in the figure. 
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Fig. 4.— Ratio method: (a) likelihoods corresponding to 201 clusters with Lx > 2.5 x lO'^^erg/s, in (b) to 397 clusters with M 500 ^ 
0.2 X IO^^Mq and in (c) to the full sample of 481 clusters. Dotted, dash-dotted and dashed lines correspond to clusters in the interval 
2 : = [0.05,0.1], 2 = [0.15,0.2] and [0.25,0.3], respectively. The red solid lines correspond to the full likelihood. Anisotropies were taken as 
averages on discs of size ^500 • 
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Fig. 5.— Fit method: Panels represent the same likelihoods than in Fig. [4| and lines follow the same convention. 



Fig. 6 .— Frequency dependence of the TSZ temperature anisotropies averaged over the cluster in the redshift subsamples. Solid lines 
correspond to the standard adiabatic evolution, a = 0, without relativistic corrections while dashed lines include these corrections. 
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Fig. 7.— Constraints on the adiabatic evolution of the background temperature at different redshifts given in Fig. [6] The results of 
Table a re represented by so lid circles with their associated error bars. For comparison, open squares and error bars correspond to the 
results of IHurier et al.l II2014I) . Points are displayed at the mean redshift of the clusters in the bin. The shaded areas correspond to the 
68 % and 99.7% confidence level. 

























